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Using the data recorded by the ANTARES neutrino telescope during 2007 and 2008, a search for 
high energy neutrinos coming from the direction of the Sun has been performed. The neutrino 
selection criteria have been chosen so as to maximize the rejection of the atmospheric background 
with respect to possible signals produced by the self-annihilation of weakly interactive massive 
particles accumulated in the centre of the Sun. After data unblinding, the number of neutrinos 
observed was found to be compatible with background expectations. The results obtained were 
compared to the fluxes predicted by the Constrained Minimal Supersymmetric Standard Model, 
and 90% upper limits for this model were obtained. Our limits are competitive with those obtained 
by other neutrino telescopes such as IceCube and SuperKamiokande, which give ANTARES lim- 
its for the spin-dependent WIMP-proton cross-section that are more stringent than those obtained 
by direct search experiments. 
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1. Introduction 



There is compelling evidence that 83% of the matter in the Universe is in the form of non- 
baryonic and non-relativistic matter that does interact weakly with the ordinary matter, the so-called 
dark matter. Much of this evidence comes from its gravitational effects on the motion of galaxies, 
clusters of galaxies and from the large scale structure of the Universe. The existence of dark matter 
is a key component of our present standard cosmological model, and results from the study of the 
CMB anisotropics and gravitational lensing in galaxy clusters further support its existence (for a 
review of the evidence, candidates and constraints, see |[T]] and [Q]). 

One of the most favoured hypothesis is that dark matter is made of weakly interacting massive 
particles (WIMPs) that are embedded in the visible, baryonic part of the galaxies and surround them 
in the form of a halo. There are a variety of candidates for WIMPs, among which those provided 
by theories based on supersymmetry (SUSY) or universal extra dimensions (UED) attract a great 
deal of interest. In some classes of SUSY, the lightest particle is stable due to the conservation of 
R-parity that forbids its decay to standard particles, making it a good candidate as a dark matter 
WIMP 

In addition to its gravitational effects, the search for evidence of the existence of WIMPs is at 
present performed, on the one hand, by looking for the recoiling products of the elastic scattering of 
dark matter particles off normal baryonic matter in suitable detectors, the so-called direct searches, 
and on the other hand, indirectly by the observation of the final products of the possible annihilation 
of WIMPs that have accumulated in astrophysical objects. WIMPs can scatter elastically in the Sun 
or the Earth and become trapped in their gravitational potential wells, accumulating in sizeable 
numbers over the age of the solar system, therefore a very wide region in the Galaxy must have 
contributed, reducing the dependence on the detailed structure of the dark matter halo. The WIMPs 
proposed by SUSY are Majorana fermions that can self-annihilate, giving rise to standard particles. 
Neutrinos will be at the end of the decay chain of the products of the WIMP self-annihilation, they 
can escape these astrophysical objects and be detected by neutrino telescopes on the Earth. 

In this paper, the indirect search for dark matter looking for high energy neutrinos coming from 
the Sun using the data recorded by the ANTARES neutrino telescope in 2007 and 2008 is reported. 
The number of neutrinos observed is compared to the neutrino fluxes predicted by the Constrained 
Minimal Supersymmetric Standard Model (CMSSM) [^], a minimal supersymmetric extension of 
the Standard Model with supersymmetry-breaking scalar and gaugino masses constrained to be 
universal at the GUT scale. 

The layout of the paper is as follows. In Section the main features of the ANTARES neutrino 
telescope and the reconstruction algorithm used in this work are reviewed. In Section ||, the Monte 
Cai^lo simulations of the signal, expected from the investigated WIMP models, and the background, 
expected from atmospheric muons and neutrinos, are described. In Section ^, the method used to 
optimise the criteria to select the sample of neutrino events and reduce the background following 
a blind data strategy is presented. Finally, the results obtained are discussed in Section ||, where 
limits on the models investigated are imposed from the absence of a dark matter signal. 
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2. The ANTARES neutrino telescope 

ANTARES is the first undersea neutrino telescope and the largest of its kind in the Northern 
Hemisphere [Q]. It is located at 2475 m depth in the MediteiTanean Sea, 40 km offshore from 
Toulon (France) at 42° 48' N and 6° 10' E. 

The telescope consists of 12 mooring detection lines made up of 25 storeys each. The standai^d 
storey is composed of a local control module that contains the front-end and slow-control electron- 
ics and three optical modules (OMs) that house a 10-inch photomultiplier. The OMs are looking 
45° downwards in order to optimise their acceptance to upgoing light and to avoid the effect of 
sedimentation and biofouling. The length of the lines is 450 m and the horizontal distance between 
neighbouring Unes is 60-75 m. The absolute time accuracy is ensured at the millisecond level by the 
UTC time provided by the GPS system connected to the clock system of the detector. The relative 
time between the elements of the detector is achieved by calibration of the lines in the laboratory 
with short light pulses, by the use of a 25 MHz clock system and by the operation of a series of 
optical beacons distributed along the lines that emit short light pulses through the water. A relative 
timing of the order of one nanosecond is reached. Additional information on the detector can be 
found in reference [Q]. 

Neutrinos are detected through the products of their interaction with the matter inside or close 
to the detector. In the channel of neutrino observation used in this work, a high energy neutrino 
interacts in the rock below the detector producing a relativistic muon that can travel hundreds of 
meters and cross the detector or pass nearby. This muon induces Cherenkov light when travelling 
through the water, which is detected by the OMs. From the time and position information of the 
hits in the OMs, the direction of the muon - which at high energy is essentially that of the neutrino 
- can be reconstructed. 

Data-taking started in 2007 when only 5 lines of the detector where installed. The full detector 
was connected in May 2008 and has been operating ever since, except for some periods in which 
repair and maintenance operations have taken place. For this work the data taken during 2007 
and 2008 has been used. Results of other searches using this data-taking period can be found 
elsewhere 

The reconstruction of the track from the position and time of the hits of the Cherenkov photons 
in the OMs is a key ingredient of the physics analysis. In this work, a dedicated fast algorithm to 
reconstruct the muon track is used [^. The algorithm is based on the minimization of a ;^^-like 
quality parameter, Q, that uses the differences of the expected and actual times of the detected 
photons coiTccted by the effect of light absorption in water. 

Monte Carlo simulations indicate that selecting tracks with a quality parameter per degree of 
freedom of the fit smaller than 1.4 results in a purity of 90% for upwai^d reconstructed multi-line 
atmospheric neutrinos, the contamination being misreconstructed downgoing atmospheric muon 
events, with an angular resolution of about few degrees at energies of tens of GeV, driven at ~ 40% 
by the kinematic of the neutrinos in low energy regime. 

3. Signal and background simulation 

The number of muon neutrinos as a function of their energy arriving at the Earth's surface 
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from the Sun's core, dNy/dEy, is computed using the software package WimpSim [Q]. The usual 
self-annihilation channels {qq, 11, WW, ZZ, Higgs doublets 00* and vv) were simulated for 17 
different WIMP masses from 10 GeV to 10 TeV. Oscillations between the three neutrino flavours 
both in the Sun and during their flight to Earth as well as V absorption and T lepton regeneration in 
the Sun are taken into account. 

For the CMSSM, three main self-annihilation channels were chosen for the lightest neutralino, 
Xi, namely: a soft neutrino channel, XiXi — ^ t>b, and two hard neutrino channels, XiXi ^ W+W^ 
and XiXi ~^ T+T^. Since which of these three channels is dominant depends on the region of the 



CMSSM parameter space being analysed [10], a 100% branching ratio was assumed for all of them 



in order to explore them on an equal footing. 

The main backgrounds for this search are muons and neutrinos produced in the interaction 
of cosmic rays with the atmosphere. Downgoing atmospheric muons dominate the trigger rate, 
which ranges from 3 to 10 Hz depending on the exact trigger conditions. They are simulated using 
Corsika [|n|]. Upgoing atmospheric neutrinos, which are recorded at a rate of ~1 mHz (about four 



per day) [^], are simulated according to the parametrisation of the atmospheric flux from [ |12| ] 
in the energy range from 10 GeV to 10 PeV. 

The Cherenkov light produced in the vicinity of the detector is propagated taking into ac- 
count light absorption and scattering in sea water. The angular acceptance, quantum efficiency and 



other characteristics of the PMTs are taken from [13] and the overall geometry corresponded to 
the layout of the ANTARES detector [Q] according to the data taking periods (from 5 to 12 line 
configurations). 

4. Optimisation of the event selection criteria 

The data set used in this analysis comprises a total of 2, 693 runs recorded between the 2T^ 
of January 2007 and the 31*' of December 2008, corresponding to a total livetime of ~294.6 days. 
The detector consisted of 5 Unes for most of 2007 and of 9, 10 and 12 lines for 2008. 

Only upgoing events are kept in the analysis. The muon tracks are required to have cos 6 < 
0.9998 in order to exclude those for which the fit stopped at the boundary. The fit is required to 
use a number of hits greater than five in at least two lines in order to ensure a non-degenerate 
5-parameter fit with a proper reconstruction of the azimuth angle. 

The UTC time of the events is uniformly randomised on the period of the data taking in order 
to estimate the background in the Sun's direction from the data itself. The local coordinates (6, 
0) are kept so as to preserve the detector geometry in the optimisation of the selection criteria. 
This procedure provides a means to follow a data blinding strategy while using all the relevant 
information on the detector's performance. 

The values used in the event selection criteria for the quality parameter, Q, and for the aperture 
of the search cone in the Sun's direction, *F, are optimised following the model rejection factor 
(MRF) technique [|l4|]. For each WIMP mass and each annihilation channel, the values of Q and *F 
used ai^e those that optimised the average upper limit on the + flux, 0v^+v^ , as defined by: 

L(^e//(,-'WwiMPj X ^eff 
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where the index / denotes the different periods in the detector infrastructure (5, 9, 10 and 12 de- 
tection lines), is the average upper limit at 90% confidence level (CL) computed from the 



time-scrambled data set and using a Poisson distribution in the Feldman-Cousins approach [|15|]; 
r^yy is the livetime for each detector configuration in 2007-2008, namely: ~134.6, ~38, '^39 and 
~83 days for 5, 9, 10 and 12 lines respectively. The effective area averaged over the neutrino 
energy, Ag//(MwiMp)', is defined as (the index / is implied): 

v,v Y Jo dE^"-'^v + dE^aHy J 

where £'(/' = 10 GeV is the energy threshold for neutrino detection in ANTARES, Mwimp is the 
WIMP mass, dNyy/dEy^y is the energy spectrum of the neutrinos or the anti-neutrinos at the 
surface of the Earth for the three channels of interest in this analysis, ?inAAeff{Eyy) is the effective 
area of ANTARES as a function of the neutrino or anti-neutrino energy. Due to their different cross- 
sections these two effective areas are slightly different and therefore ai^e studied separately. 



5. Results and discussion 

Once the optimised values of Q and *F were obtained using the time-scrambled data, the data 
sample is unblinded. Figure |I] shows the distribution of the spatial angle between the tracks of the 
events and the Sun's direction obtained after applying the basic selection criteria on the zenith angle 
and the minimum number of hits and lines. A total of 27 events were found within a 20 degrees 
spatial angle and no excess in the Sun's direction above the scrambled background are observed. 

Using the values for the cuts obtained in the optimisation procedure, limits on the + 
flux, 0v^+v^» can be extracted from the data according to: 




» 10^ 10= 10 

Mw,„p(GeV) 



Figure 1: Left: Distribution of the spatial angle m g[0",20''] of the event tracks with respect to the Sun's 
diretion for the expected background computed from the time-scrambled data (solid blue line) compared 
to the data after the basic selection criteria (black triangles) .Ala Poisson uncertainty is shown for each 
data point (black cross). Right: 90% CL upper limit on the muon flux as a function of the WIMP mass 
in the range Mwimp e[50 GeV;10 TeV] for the three channels bb (green), VK+W" (blue) and T+T^ (red, 
ANTARES only). The results from SuperKamiokande 1996 - 2008 [[7| (dotted lines) and IceCube-40 plus 
AMANDA 2001-2008 (dashed lines) are also shown. 
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^90% 



where ji^^'^'' is the upper Umit at 90% CL on the number of observed events and the rest of variables 



have the same meaning as in Eq. 4.1 



The conesponding limits for muons are calculated using a conversion factor between the neu- 
trino and the muon fluxes (^^ = Fy^p x 0v^+v^) computed using Dai^kSusy [16]. Figure [l| (on the 
right) shows the 90% CL muon flux limits 0^ for the channels bb, W^W^ and t+t^. The latest 
results from Super Kamiokande [ p^ ] and IceCube-40 plus AMANDA [18] are also shown for com- 
parison. Despite its smaller detector volume, the ANTARES limits both in the soft and hard chan- 
nels are similar to those obtained by IceCube40 and AMANDA in the mass range Mwimp €[50; 100] 
GeV. In this mass range most of the sensitivity of the South Pole telescopes comes from AMANDA- 
II, whose effective area and energy threshold were similar to those of ANTARES in its 12-line 
configuration. Compared to the limits set by SuperKamiokande, those of ANTARES are more 
stringent in the high mass region Mwimp > 150 GeV. 

Systematic uncertainties have been taken into account and included in the evaluation of the 
limits using the approach of reference [ [l9| ] by means of the Pole software. The total systematic 
uncertainty on the detector efficiency is around 20% and comes mainly from the efficiency and time 
resolution uncertainties of the OMs, the total angular resolution and the absolute pointing accuracy. 
This uncertainty translates into an increase of the upper limit between 3% and 6% depending on 
the WIMP mass. 



SI 




Mw™p(GeV) Mw,„p(GeV) 



Figure 2: 90% CL upper limits on the SD and SI WIMP-proton cross-sections (left and right hand respec- 
tively) as a function of the WIMP mass in the range Mwimp €[50 GeV; 10 TeV], for the three channels: bb 
(green), VK+W^ (blue) and t+t^ (red), for ANTARES (solid line) compared to the results of other indi- 
rect search experiments: SuperKamiokande 1996 — 2008 ^% (dotted lines) ,IceCube-40 plus AMANDA 
2001-2008 [111 (dashed lines) and Baksan 1978-2009 [g (colored dot-dashed line) and the result of 
the most stringent direct search experiments (black): PICASSO 2007 — 2012 | ]23| ] (three-dots-dashed line), 
KIMS 20 11 - 2012 [01 (two-dots-dashed line), SIMPLE 2004 - 201 1 [Q (one-dot-dashed line), COUPP 
2010 - 201 1 [Hi (one-dot-large-dashed line) and XENONIOO 201 1 - 2012 ^ (dashed hne). For a com- 
parison to the theoretical expectations, grid scans corresponding to the CMSSM model have been added. 



Assuming equilibrium between the WIMP capture and self-annihilation rates in the Sun, the 
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limits on the spin-dependent (SD) and the spin-independent (SI) WIMP-proton scattering cross- 
sections can be obtained for the case in which one or the other is dominant. 

A conservative approach to the dark matter local halo was considered assuming a gravitational 
effect of Jupiter on the Sun's capture rate, which reduces about 13%(87%) the SD capture and 
about 1% (12%) the SI capture for 1(10) TeV WIMPs. A local dark matter density of about 0.3 
GeV.cm^^ was assumed. No additional dark matter disk that could enhance the local dark matter 



density was considered (see [20] for a discussion). 



The 90% CL limits for the SD and SI WIMP-proton cross-sections extracted from the chan- 
nels bb, W^W^ and T+T^ are presented in Figure ^ The latest results from SuperKamiokande, 
IceCube-40-i- AMANDA and Baksan together with the latest and the most stringent limits from the 
direct seaixh experiments PICASSO, KIMS, SIMPLE, COUPP and XENONIOO are also shown. 
The allowed parameter space from CMSSM model according to the results from an adaptative grid 
scan performed with DarkSUSY [ [l^ ] are also shown. For a proper comparison all the limits pre- 
sented in Figure § are computed with a muon energy threshold at f"^ = 1 GeV. For these figures, 
the shaded regions show a grid scan of the model parameter space taking into account the last limits 
for the Higgs boson mass from ATLAS and CMS merged together such as Mh = 125 ±2GeV [||]. 
A relatively large constraint on the neutralino relic density < Q.cDMh^ < 0.1232 is used not to 
enclose the studied dark matter particle to only one possible nature. 



The neutrino flux due to WIMP annihilation in the Sun is highly dependent on the capture 
rate of WIMPs in the core of the Sun, which in turn is dominated by the SD WIMP-proton cross- 
section. This makes these indirect searches better compared to direct search experiments such as 
KIMS and COUPP. This is not the case for the SI WIMP-proton cross-section, where the limits 
coming from direct search experiments such as CDMS and XENONIOO are better thanks to their 
target materials. Therefore, there is a sort of complementarity between both types of searches. 
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